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ABSTRACT 



We study irreversible A-B reaction kinetics at a fixed interface separating two immiscible 
bulk phases, A and B. Coupled equations are derived for the hierarchy of many-body corre- 
lation functions. Postulating physically motivated bounds, closed equations result without 
the need for ad hoc decoupling approximations. We consider general dynamical exponent 
z, where x t ~ t x l z is the rms diffusion distance after time t. At short times the number 
of reactions per unit area, R t , is 2nd order in the far-field reactant densities n'^,n'§. For 
spatial dimensions d above a critical value d c — z — 1, simple mean field (MF) kinetics per- 
tain, R t ~ Qbtn'^n'^ where Qb is the local reactivity. For low dimensions d < d c , this MF 
regime is followed by 2nd order diffusion controlled (DC) kinetics, Rt ~ xf +1 n°£n^ , provided 
Qb > Qt ~ (rig > )[ z_ ( d+1 )l/ d . Logarithmic corrections arise in marginal cases. At long times, 
a cross-over to 1st order DC kinetics occurs: R t ~ x t n A . A density depletion hole grows on 
the more dilute A side. In the symmetric case (n^ = rig), when d < d c the long time decay 
of the interfacial reactant density, n s A , is determined by fluctuations in the initial reactant 
distribution, giving n s A ~ t^ d /( 2z ). Correspondingly, A-rich and B-rich regions develop at 
the interface analogously to the segregation effects established by other authors for the bulk 
reaction A + B — > 0. For d > d c fluctuations are unimportant: local mean field theory 
applies at the interface (joint density distribution approximating the product of A and B 
densities) and n s A ~ t( 1 ~ z )/( 2z ). We apply our results to simple molecules (Fickian diffusion, 
z = 2) and to several models of short-time polymer diffusion [z > 2). 



PACS numbers: 05.40. +j, 68.45.Da, 82.35.+t 
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1. Introduction 



In a large class of chemically reacting systems, irreversible bimolecular reactions occur at a 
permanent interface separating two bulk phases. Reactive molecules in one phase are able to 
react with molecules in the other phase only; hence reaction events can occur within the limits 
of the narrow interfacial region only. A number of technologically important examples |T]fJ] 
entail small molecules reacting at liquid-liquid, liquid-gas or liquid-solid interfaces. In the 
present study we address interfaces which are fixed in space and do not broaden as reactions 
proceed; the two bulk phases do not mix with one another. However, the physics we will 
explore may provide insight to the very different problem of non- stationary reactive fronts 
where chemical reactions occur at a moving and possibly broadening interface separating 
miscible phases @,|g§rj,|,0,|IUig]. 

Another important class of reactive interfacial systems involves macromolecules. Of par- 
ticular technological significance is the process of reactive blending [[L^IH] where the compat- 
ibilization of two immiscible polymer melts A and B is assisted by attaching reactive groups 
to a certain fraction of the chains. The A-B copolymers generated by reactions, which can 
occur at the A-B melt interface only, serve both to reinforce the interface ||bl| , |l5f and to 
promote the mechanical mixing of the two melts |T3| , |T6| . 

The manner in which reaction kinetics are modified by the presence of an interface is a 
fundamental issue within the general field of chemical reaction kinetics. Despite this, and 
despite the numerous applications such as those mentioned above, no complete and system- 
atic theory exists. We emphasize that each reaction event necessitates the simultaneous 
arrival, at the same location within the interface, of two molecular species A and B, one 
from each bulk phase. (This is very different to the problem [|T7| , |r8f of a single bulk adjacent 
to a homogeneous "reactive interface" where each "reaction" event, e.g. the irreversible ad- 
sorption of a molecule onto a solid surface, requires the arrival of only one molecule at the 
interface.) The interfacial reaction kinetics which are the subject of the present paper were 
theoretically studied for the case of small molecules by Durning and O'Shaughnessy and 
the end-functionalized polymer case by O'Shaughnessy and Sawhney and Fredrickson 

p2| . These theories in fact described a certain short time regime only, for systems where the 



reacting species are dilute in an unreactive background. Fredrickson and Milner |23j argued 
that at later times different kinetic behaviors onset, with forms dependent on reactive species 
concentration. 

In this paper we develop a near-exact theory of irreversible interfacial reaction kinetics. 
We calculate time-dependent reaction rates as a function of density and local reactivity of the 
reactive species. In addition, density profiles on either side of the interface are determined. A 
short version of the present manuscript has appeared ||24|| . Our framework is quite general in 



terms of the diffusive dynamics of the reactive species, as defined by the dynamical exponent, 
z : 

x t = a y , (1) 

where x t is the rms displacement of a reactive group after time t. Here a is the linear 
dimension of the reactive species A and B, and t a is the diffusion time corresponding to a. 
Thus, setting z = 2 in our results yields reaction rates for small molecules obeying Fick's 
diffusion law. As a second example, if one seeks the short time reaction kinetics of small 
reactive groups attached to polymer chains in the melt, then appropriate values would be 



or z = 8, depending on the time regime and degree of entanglement |25| , |26 
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This study will always assume the interface is clean. Thus we ignore effects associated 
with accumulation of A-B reaction products at the interface whose presence may eventually 
diminish reaction rates [E0lETI,E3 . 



Bulk Kinetics (A+B — > 0): a Brief Review 

Before attacking the present interfacial problem, it is helpful to consider first the analogous 
and somewhat simpler problem in the bulk, where irreversible bimolecular reactions between 
A and B, generating inert products, can occur anywhere within a single bulk phase. There 
are no A- A or B-B reactions. For the small molecules case, z = 2, many well-established 
results exist ( U A + B — > 0") [p^, p^ , p^j30|j3T|1 . Suppose A and B have equal diffusivities and 



are initially uniformly distributed with equal densities n(0) within some solvent. (The case 
n(0) = l/a d would correspond to every molecule being reactive; generally n(0) < l/a d .) 
Reactions are now "switched on" at t — 0. Then, whenever an A and a B particle collide 
(i.e. approach to within distance of order a of one another) they react irreversibly with 
probability Qb per unit time, where Qb is the local reactivity. For simplicity, let us confine 
our bulk discussion to "infinitely" reactive particles; reaction then occurs every time an A- 
B pair collides. This corresponds to setting Q b = l/t a (the effective local reactivity 
cannot exceed the rate, l/t a , at which diffusion can bring two reactive species together). 
Our discussion will consider a general spatial dimensionality d. 

What are the reaction kinetics in this bulk system? What are the time dependencies of 
the number of reactions per unit volume which have occurred after time t, namely TZ^ nlk , 
and the reactant density n(t)? The simplest guess is that mean field (MF) theory applies: 
this amounts to assuming reactants are always distributed as in equilibrium. Hence the 
reaction rate equals the equilibrium density of A-B pairs in contact, multiplied by Qb- Thus, 

7t 4 bulk = (d/dt)n]; ulk = Q b a d n 2 {t) = (a d /t a ) n 2 {t). Now this MF prediction is in fact valid 
only for sufficiently large d such that diffusion is effective in dissipating reaction-induced 
non-equilibrium spatial correlations. The maximum number of A-B pairs which diffusion 
can have brought together by the time t increases as xf ~ t d l z \ provided d > z, this is 
sufficient to restore the depletion, which would arise in the two-body correlation function, 
due to reactions as implied by the MF prediction, 7^ ulk ~ t. But for lower dimensions, 
d < z, since diffusion cannot supply material fast enough to keep pace with this reaction 
rate, equilibrium spatial correlations are destroyed: a depletion hole of size x t grows in the 
A-B 2-body correlation function. Reaction kinetics are then very different; for short times 
-^buik ^ x d n 2(Q^ j s the number of reactive pairs initially within diffusive range of one another, 
i. e. whose initial separations were less than xt- In summary, for times short enough that the 
relative density drop is small, n{t) w n(0), we have 



a d /t a (d > z) 

dxf/dt ~ t d l z ~ l (d < z) 



7t t bu * = = k hulk n 2 (t) , A; bulk (t) « { (short times). 

dt 



(2) 

These are second order rate kinetics, with a rate constant, k hnlk , which is time-dependent 
for low dimensions d < z. 

The two classes of kinetics in eq. (|2|) reflect the fact that reactive groups explore space 
"compactly" in low dimensions |33|j3^ , |35|| : for d < z it is simple to show that (in the 
absence of reactions) the number of collisions between an A-B pair, with some given initial 
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separation, increases for large times as ~ t x ~ d l z . Thus reaction is inevitable by the time t for 
any reactive pair initially separated by x t or less. By contrast, for d > z space is explored 
in a "non-compact" or dilute fashion; with finite probability, the same two particles may 
avoid collision as t — > oo. This survival probability is an increasing function of the initial 
separation. Reaction is no longer inevitable between all pairs within diffusive range, and 



MF theory applies [p6| , p2| . For the interface problem, we will find a similar division between 
compact and non-compact reaction kinetics, but now at a dimension d + 1 = z. Indeed, the 
short time interface kinetics turn out to be analogous to those of a d + 1-dimensional bulk 



problem PJ20J21 



Eq. (Q) describes the short time kinetics. What happens at very long times? For d > z 
one might anticipate the MF kinetics of eq. (0) would continue indefinitely, implying n ~ 1/t 
asymptotically. In fact, for two-species A-B systems this is true only for very high dimensions, 



d > 2z. This was demonstrated for z = 2 by Ovchinnikov and Zeldovich [271 and by 



Toussaint and Wilczek pSj . These authors showed that in lower dimensions fluctuations 
in the initial density distribution determine the asymptotic form of n(t). Their argument, 
generalized to arbitrary z, is roughly as follows. Consider a portion of the reacting system 
of volume Q and let iV^(t), Ns{t) be the number of unreacted A and B particles in this 
region at time t. Assuming random initial spatial distributions of A and B, the initial 



fluctuations of A^(0) and N B (0), about their mean value n(0)Q, will be of order yn(0)Q. 
(Note that n(t) is the mean density after time t.) Of the same order will be the fluctuations in 
8 No = Na{0) — iVg(O), the average value of which is zero. As reactions proceed, fluctuations 
will diminish. However, since reaction events conserve the difference between the number of 
A and B particles, fluctuations in 8N t = Nj\,(t) — Nsit) can decay by diffusion only. Thus, if 
we consider small regions, Q < xf , then by time t the initial difference of order 8Nq has had 
sufficient time to decay away due to diffusion. But for large regions, Q > xf, the difference 
must be close to its original value, 8N t « 8N . Roughly, then, in a region of volume xf, 

the total number of reactants cannot be smaller than a number of the order of \Jn(0)xf. 



It follows that the density cannot decay faster than n{t) « ^Jn(0)xt/xf ~ Jn(0)t d /( 2z \ 
For d < 2z, this is a slower decay than the MF t~ l prediction, and one concludes that this 
diffusive relaxation of initial fluctuations then governs the long time decay. To summarize, 

f r 1 (d > 2z) 
n(t -> oo) ~ I (bulk). (3) 

[ r d / 2z (d < 2z) 

For interfacial reactions, we will establish a rather similar long time fluctuation-dominated 
decay of densities near the interface for sufficiently small d. Analogously to the bulk case, 
this is accompanied by segregation of reactants into A-rich and B-rich domains of size x t in 
the region adjacent to the interface. 



Interfacial Kinetics: Scaling Arguments 

Let us turn now to the interface problem, shown schematically in fig. [I]. We consider 
two d-dimensional bulk phases containing, respectively, reactive species A and B with initial 
densities n°£ and n^. The reactants are of size a < h, where h is the width of the thin (d- 
l)-dimensional interfacial region which is the locus of all reaction events. We assume A and 
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B have identical diffusion dynamics. To begin, consider the symmetric case n°£ = = n 
and the infinitely reactive limit, Qb — > l/t a (every A-B collision produces a reaction). 

The short time reaction kinetics of this interfacial system are analogous to those of a 
(d + l)-dimensional bulk problem. To see this, consider how many degrees of freedom are 
needed to specify the "reaction rate" for a single A-B pair. One coordinate must specify 
how far from the interface particle A lies, and similarly for B. A further d — 1 coordinates 



must specify their relative location, giving d + 1 degrees of freedom [19] in total. That is, 
there are d + 1 diffusive degrees of freedom which must vanish in order that an A-B pair 
may react. These are the reaction conditions for a d + 1-dimensional bulk diffusion-reaction 
problem, and similar reasoning to that for the bulk dictates that non-compact MF kinetics 
pertain for d + 1 > z, whilst for d + 1 < z the kinetics are of compact diffusion-controlled 
(DC) form. Thus the reaction rate per unit interface area, TZt, obeys 2nd order rate kinetics 
with a 2nd order rate constant k^ : 

( ha d /t a (d + 1 > z) 

■JZt = k (2) n 2 , k {2 \t) « I (short times , Q b t a = 1) 

( dx d+1 /dt ~ tV+W*- 1 (d + Kz) 

(4) 

in complete analogy to eq.(0) for the bulk problem, but with d replaced by d + 1. The mean 
field result for d+1 > z follows because in equilibrium there are ha d n 2 A-B pairs in contact 
per unit area of interface The DC compact kinetics are determined by the small fraction 
(at short times) of A-B pairs which were initially separated by less than x t ; for d + 1 < z 
any such pair will definitely have reacted by time t. The number of such A-B pairs per 
unit interfacial area is xf +1 n 2 (see fig. p. Note that the dimensions of k^ are x d+1 /t, as 
appropriate to a {d + l)-dimensional bulk problem. For the remainder of this paper we will 
refer to d + 1 > z and d + 1 < z as the "non-compact" and "compact" cases, respectively. 

Consider the long time behavior now. This is completely different to the bulk. In fact at 
long times the effective dimensionality of the problem changes from d+ 1 to 1 and, moreover, 
the reaction kinetics become of first order. Let us first investigate the compact case, d+1 < z. 
Consider an A particle that was initially within a distance I of the interface, where I = n~ x l d 
is the typical separation between reactants. By time ti = t a (na d )~ z ^ d , its exploration volume 
will typically contain one B particle, in the other bulk, which was initially within / of the A 
particle. Since d + 1 < z, reaction is certain. It follows that for times t > ti the interface 
becomes, in effect, "perfectly absorbing:" almost every reactive species reaching it will suffer 
a reaction. Thus a density depletion hole develops at the interface (see fig. ||) and the 
reaction rate is limited by diffusion to the interface, TZt ~ %tn. One concludes that long time 
reaction kinetics are now first order, with a first order rate constant k^ given by 

TZ t = k^n, k^^^r^t 1 ^ 1 (*-oo) . (5) 

Notice that the dimensions of k^ are x/t, as would be appropriate to a one- dimensional bulk 
problem. In correspondence to the kinetics being first order, this DC regime is accompanied 
by a growing hole of size x% in the one-body density "correlation function," i. e. in the density 
field itself, n(r). This is very different to the hole, also of size x t , which grew in the two-body 
density correlation function for the second order d + 1 < z compact DC kinetics at short 
times, eq. (|j). For that regime, the density field itself was unchanged from equilibrium. 

1 A slight complication here, which will be addressed in section 6, is that this result is modified when 2<d+l<z + l;in 
that case fc( 2 ' h d+1 /tf l where = t a (h/a) z . 
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What are the long time kinetics for the non-compact case, d + 1 > zl The answer is: 
the same as for the compact case. The only difference is that the cross over from d + 1 to 
1-dimensional behavior no longer occurs at ti, but at a timescale we name t* m . In this case 
we can estimate t* m using a mean field picture since the early kinetics are MF. Consider an A 
particle initially within xt of the interface, as in fig. |2|. After time t, it has made t/t a "steps," 
a fraction (h/xt) of which were within the interfacial region where B particles are present at 
density n. Thus for each of these interfacial steps the probability that the A particle was in 
contact with any B particle is na d . Hence the total reaction probability, P m , is 

P m (t) « (-f) (na d ) « - na d (t/t.) 1 - 1 '* (Q b = l/t a ) . (6) 

Setting P m {t*m) = 1, one obtains t* m /t a = [a/(hna d )) z ^ z ' 1 \ Thus, for t > t* m any A particle 
within diffusional range of the interface will definitely have reacted with a B. This is a many- 
body effect; by t* m any A near the interface is bound to have reacted due to the mean reaction 
field created by all of the B molecules. We conclude that for large times a density depletion 
hole develops also for the non-compact case, following the same kinetics as eq. (^j). 

So far we have considered "infinitely" reactive species, Qb ~ l/t a - Such local reactivities 
Qb are realized for radicals f37jl , and in certain other processes such as phosphorescence 



quenching However, these are very exotic exceptions to the general rule: in virtually all 
practical situations Qb is tiny, Qbt a ~ 10~ 6 . Indeed, for the vast majority of reacting species, 
Qb values are many orders of magnitude smaller than |39|.|4(|. It is essential, therefore, 
to establish how the picture we have developed above is modified for finitely reactive systems. 
For non-compact cases, d + 1 > z, there is no qualitative change from the kinetics of eq. ([|): 
again, a short time 2nd order MF regime, now with = Qbha d , is followed at t* m by a DC 
first order regime, but now the formula for t* m is modified. Notice that the expression for 
P m (t) in eq. (^) is the mean number of collisions experienced by the A particle; multiplying 
this by the reaction probability per collision, Qbt a , yields the total reaction probability for 
general Qb- Defining P m (t£j = 1, we have 

, / 1 \ 1 -^ z t* ( 1 \ z/( - z - 1) h 

P m (t) w (Qt a ) na d ( — J , J± = U—\ , Q = Q b (7) 

where Q emerges as an effective local reaction rate coarse-grained over the interface width 
h. 

In the compact case, d + 1 < z, kinetics are more fundamentally modified by finite 
reactivity. For Qbt a = 1 we have seen an initial 2nd order DC regime followed at f by 1st 
order DC kinetics. But for Qbt a < 1, a new MF regime appears at early times. Consider an 
A-B pair near the interface whose members are initially closer than xt to one another, as in 
fig. |l|. By time t, A has taken (t/t a )(h/xt) steps in the interface. For a fraction (a/xt) d of 
these, B was in contact with A since B is equally likely to be anywhere within its exploration 
volume xf. Thus, the 2-body reaction probability for this pair is given by 

This implies a characteristic timescale, t%, defined such that ^2(^2) = 1) 



P2(t)*l-)[z3)(r)ma)*Qt a ( r ) • (8) 
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For t > t\ any pair initially within diffusive range will definitely have reacted; this tells us 
that kinetics must have 2nd order DC form for such times. Thus the DC regime of eq. (H) 
begins only at t\\ for shorter times, t < t%, since -P^) "C 1, correlations are little disturbed 
from equilibrium and it follows that MF kinetics apply, = Qbha d . In fact, for sufficiently 
small Q ("weak systems") t\ will exceed t* m in which case the 2nd order DC regime will 
disappear. In later sections we will carefully distinguish between this case and the case of 
"strong systems" (t* 2 <t* m )- 

Interfacial Kinetics: the Technical Difficulties 

In the present work we will develop a near-exact formalism to justify these scaling argu- 
ments. The difficulty is the many-body character of this problem. Consider for example the 
reaction rate per unit area, 1Z t . This equals the number of reactive A-B pairs per unit area 
which are in contact at the interface, p^g(i), multiplied by the local reactivity Qb'- 

^ = Xp s AB (t) , A = Q b ha d = Qa d+1 (10) 

where the quantity A will turn out to be a natural coupling constant in our theory. Now 
p\ B (t) is the two-body density correlation function pab(*a, r B', t) (the number of A-B pairs 
at ya,yb per unit volume squared) evaluated at the interface, ya = yb = 0: 

P S A B (t) = PA B (0,0;t) . (11) 

We take the origin of our coordinate system to lie on the interface plane and we have used 
translational invariance in the directions parallel to the interface (hence Pab(^) * s spatially 
uniform). One sees that to determine the reaction rate we need information on the two- 
body density correlation function. However, any dynamical equation for the latter inevitably 
involves three-body correlation functions Paba, Pbab- The dynamics of these objects in turn 
involve higher order correlations, and so forth. This hierarchical structure is the signature 
of the many-body nature of the problem. 

How can a theory deal with these many-body complexities? One possible approach |23[ , 
a mean field approximation, would be to assume p% B (t) = [n s (t)} 2 , where n s (t) ee n(r = 0) is 
the density of A (or B) reactants at the interface. This approximation cannot always be valid: 
for example, in the compact case, d+l < z, this would disagree with the short time 2nd order 
DC behavior of eq. since in this regime the density field is unchanged from equilibrium, 
n s (t) ~ n(0); hence the assumption p s AB {t) = [n s (t)} 2 would wrongly yield 1Z t ~ t. Does 
this approximation make sense at longer times? Now since we have established (eq. (j5|)) 
the asymptotic law TZt ~ nt^~ z ^ z , this approximation would then imply n s (t) ~ t( 1 ^ z )/( 2z ) 
which as we will see is correct for the non-compact case only. For the compact case, the 
long time decay of n s (t) is in fact determined by the rate at which fluctuations in the initial 
distribution of A and B reactants decay. This gives rise to a different decay law, invalidating 
the local mean field approximation. 

To see how densities at the interface, n s (t), decay for large times, consider a simple 
generalization of the argument of Ovchinnikov and Zeldovich and Toussaint and Wilczek, 
extended to the interface problem. Consider a region at the interface of volume Q, half of 
which is on the A side and half on the B side. The difference 5N(t) = N A {t) — N B (t) between 
the number of A and B in is initially of order vn?2. Now fluctuations in 5N(t) can decay 
by diffusion only. Only if Q is smaller than xf did these fluctuations have sufficient time 
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to have decayed by time t. For bigger regions, 8N(t) ~ 5N(0) « y/n£l. Thus reactant 



densities at the interface, for example, cannot decay faster than Jnxf/xf ~ \fnt d /( 2z \ In 

the compact case, d + 1 < z, this is a slower decay than [Pab^)] 1 ^ 2 ■ Thus the local mean 
field assumption is wrong, and subtle correlations between reactants determine the long 
time decay. Correspondingly, for the compact case only, there is a segregation of reactants 
adjacent to the interface into A-rich and B-rich regions of size x t . 

Various approximation schemes have been used to treat reaction kinetics in the bulk. 
Typically, the three-body density correlation function is truncated in terms of lower order 
correlations; this reduces the hierarchy of reaction-diffusion equations for the many-body 
correlation functions to a closed set which are solved numerically (see ref.ffT] and references 
therein). The ad- hoc nature of such approximations is balanced by their success, as judged 
from direct numerical simulations [H|. Rigorous analysis was initiated by Doi |52| who 



developed a general formalism mapping classical many-particle systems onto quantum field 
theoretic models. Doi's formalism has been the starting point of recent renormalization 
group approaches to bulk reacting systems ||3|,|44|,|31 . 



Our approach is rather different to previous ones. We make a small number of simple 
assumptions which on physical grounds we believe are correct: we assume bounds on certain 
density correlation functions, and we assume the reaction rate to be a decreasing function of 
time which is asymptotically a power law. It is possible that these bounds might be proved 
rigorously, but we do not attempt this here. Having made these assumptions, the subsequent 
analysis is exact. In the case of systems such as reacting polymers which are not point-like 
(all the internal polymer degrees of freedom are involved in addition to the locations of the 
reactive groups) and for which z 7^ 2 at small times, our analysis, though not exact, provides 
a framework for calculating all physically interesting quantities. 



The rest of this paper aims to justify the scaling arguments presented above. In Section 
2 we present an exact mathematical formulation of the problem. We significantly simplify 
the problem in section 3 by postulating bounds on a three-body density correlation function. 
This allows us to solve for the reaction rate. In Sections 4, 5 and 6 we solve for the reaction 
rate in the compact, noncompact and marginal (z — d + 1) cases. Our results verify the 
scaling arguments presented above. In Sections 7 and 8 the density profile is calculated, 
including fluctuation effects and reactant segregation. We conclude with a discussion of our 
results in Section 9. 



2. Interfacial Pair Density, p s AB 

According to eq. (|I0|) the reaction rate is proportional to the density of A-B pairs which 
are in contact at the interface, Pab^)- I n this section, we will obtain an exact self-consistent 
integral expression for p s AB . 

We consider the general situation, illustrated in fig. [I], where the initial reactant densities 
n< Ai nr B are n °t necessarily equal. (The entire discussion of section I treated the symmetric 
case = for simplicity.) Our convention will always be that > n^. We choose the 
d-dimensional A and B bulk phases to occupy x > and x < respectively, with x being 
the direction orthogonal to the interface, and we assume that A and B species have identical 
dynamics (i.e. dynamical exponent z). 
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Throughout this paper, we use the convention that superscript T denotes a d- dimensional 
vector lying in the (d-l)-dimensional interface. Thus by definition the x-component of r T 
vanishes. 

Let us begin by treating the case z = 2, which is then simply generalized to arbitrary z. 
The second-quantization representation for classical many-particle systems developed by Doi 
|42| and by Zeldovich and Ovchinnikov ||5| allows us to derive an exact reaction diffusion 



equation for the two-body correlation function Pab(?Ai r B', t). Using Doi's formalism, we 
show in Appendix A that for small non-interacting Fickian molecules (z = 2) with diffusivity 
D 

j - D [V 2 A + p AB (r A ,r B ;t) = -X5(x A )5(r A - r B ) p AB (r A ,r B ;t) 

- \ 5(x A ) p ABB (r A ,v B ,r A ;t) - X5(x B ) Paba(j^a-> r B> r Bi t) , 

(12) 

with reflecting boundary conditions at x — 0. Note the appearance of the coupling constant 
A = Qbha d introduced in eq. fllO|). The 3-body correlation function p ABB (r A ,r B ,r' B ;t) is 
the probability density to find an A-B-B triplet at locations r A ,r B ,r' B . A similar definition 
applies to p ABA (r A , r B ,r' A ;t). 

The sink terms on the right hand side of eq. (12) describe the three ways in which 
reactions can diminish Pab{ y a, r B] t). (1) The first two-body sink term represents reactions 
between A-B pairs located at r A ,r B . The delta functions restrict reactions to r A ,r B values 
such that both A and B are in contact (i. e. within a of one another) and both A and B 
are within the interface of width h located at x — 0. These restrictions introduce a factor 
ha d . This is a somewhat coarse-grained description: our "minimal" delta-function sinks are 
appropriate provided we avoid timescales of order th = t a {h/a) z or smaller. (2), (3) The 
remaining two sink terms in eq. (12) describe reactions involving just one particle of an 
A-B pair at t a ,t B - Such a reaction involves a third particle, weighted by the appropriate 
3-body correlation function. These are many-body terms; were they absent, one would have 
a relatively simple closed 2-body system. In the next section we will deal with this difficulty 
by assuming bounds on the forms of these 3-body correlation functions. 

Consider a general value of z now, for which the two particle free propagator is G t (r A , r' A ] r B , r' B ), 
namely the probability density an A-B pair is at r^, at time t given initial location r^, r' B , 
in the absence of reactions. Without reactions A and B particles are statistically indepen- 
dent; thus G t can be written as a product of single particle propagators : 



G t (r A ,r' A ;r B ,r' B )=Gi 1 \r A ,r' A )G ( t 1 \r B ,r' B ) . (13) 



Since G^jp has only one characteristic scale, x t} dimensional analysis dictates the scaling 
form 

i f f(Ux,v x ) (|u- v| < 1) 

Gf (r, r') = - g{r/x t , v'/x t ) , g(u, v) - J (14) 

Xt { (|u-v| > 1) 

where f(u x ,v x ) is a function of order unity for every value of its arguments (u x and v x are 
the x components of u, v, respectively). The fact that / depends on u x ,v x is a result of the 
broken translational invariance in the x direction due to the reflecting boundary at x — 0. 
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Returning to eq. (12), we can write a self-consistent expression for Pab{ y a, r B]t) in terms 
of the free propagator G t . Setting = r B = one obtains 



p AB (t) = n%n% - X I dt'S^Xt - t>)p\ B {t>) - X[li{t) + I*{t)] , (15) 

where 



o 



S^\t) ee f dr T 'G t (0, r T '; 0, r T ) » (16) 

x t 

is the two-body "return probability," namely the probability density an A-B pair is in contact 
at time t at the interface, given it was in contact somewhere within the interface at t — 0. 
We have used eq. ([14]) to show that S( d+1 ^(t) ~ l/x d+l has the same scaling form as the 
return probability in a (d + l)-dimensional bulk problem. In eq. (|T5|), the two-body integral 
involving S^ d+1 \t) represents depletion in the interfacial reactive pair density p s AB (t) due to 
A-B pairs whose members reacted with one another at times t' < t and therefore failed to 
reach the origin at t (see fig. £|). The terms i^(t), Im(t) measure depletion due to many-body 
effects: 

li{t) = J o dt' J dr^dr' B G t ^(0,r^;0,r' B )p ABB {r^,r' B ,r^;t') , 

Cit) = j|Vy dr' A drl'G t ^(0,r' A ;Oy B )p ABA (r' A ,rl\r T B -,t') . (17) 

These integrals subtract off any A-B pair only one member of which was involved in a 
reaction an earlier time (see fig. |j). 

We will see later that the two-body integral involving ( S , ( d+1 ) (t) in eq. (|l^) is important at 
short times; for such times reaction kinetics are hence like those in a (d+ l)-dimensional bulk 
reaction problem. At longer times the many-body terms I^(t), I^(t) are always dominant 
and kinetics cross over to one-dimensional form. 

Eqs. (|15D, Ql6D and (17) are immediately generalized to arbitrary dynamics with arbitrary 
values of z: one simply replaces the Gaussian (z = 2) propagator G t , describing Fickian 
diffusion, with the appropriate propagator describing the dynamics. Now this would be a 
true statement for the abstract concept of small (i. e. point-like) molecules obeying x t ~ 
with z ^ 2. However, in practice non-Fickian diffusion normally results from the small 
reactive species belonging to a large structure with complex internal dynamics. The most 
important case is when the reactive group is a single monomer unit belonging to a polymer 
chain of N units. In these cases the dynamics of eq. ([15]) are not exact because they 
incorrectly presuppose a closed relationship in terms of the degrees of freedom specifying 
the location of the reactive species only. A proper treatment must first average out the 
other degrees of freedom (e.g. the locations of the other N — 1 monomers in the polymer 
case); this is non-trivial and requires renormalization group (RG) methods [EH]. However, 



RG studies of 2-body bulk polymer reaction kinetics p6| , |36| indicate that the basic physics is 
completely captured by the approximate closing of the system in terms of these coordinates 
only: correct scaling behaviors are obtained, only the prefactors being unreliable. These 
issues are discussed in detail in refs. [T7[[|6[ , P51 Therefore for the remainder of this paper we 
assume the validity of eqs. (f[5|), (0) and (17) for any value of z. 
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3. Structure of Many-Body Integral Terms l£(t),I%(t) 



In Section 2 we derived a self-consistent solution for the interfacial reactive pair density 



p% B (t), eq. dl5|) . Unfortunately this is not in a closed form for p% B (t), since the many- 
body terms involve higher order correlation functions. In this section we introduce our three 
simple, physically motivated assumptions. These enable us to express I^(t), I^(t) in terms 
of p\ B {t), which in the following section will allow us to obtain a closed solution for p s AB . 
Most calculational details will be left for Appendix B. 

Consider the three-body correlation function pabb appearing in I^(t) of eq. (17). Let us 
introduce the conditional density of B particles at r B given an A-B pair at the origin, 

/ m n +\ - PABs(0,r B ,0;t) 
pBAB{r B 0,0;t = — . 18 

PabW 

Noting that translational invariance parallel to the interface plane allows the replacement 
Pabb(ta> r B> v a ; f ) Pabb(0, t' b - y t a \ 0; f), we can express l£(t) as 

Iiit) = J*dt'J dv' B F t ^{v' B ) p BAB (r' B \0,0;t') p s AB (t') , 

F t _ t ,(r' B ) = J dv T A ' G t _t/(0, r T A '; 0, r' B + ) • (19) 

In fig. || we identify two physically distinct space-time regions which contribute to I^(t) 
in the r' B ,t' integration of eq. (19). The assumptions we are about to introduce are based 
on the following expectations about the behavior of the conditional 3-body density in these 
two regions. In region I, defined by points with x-coordinate x' B > x t ', the conditional 
density at time t' approximates its far field value, PBAB( r ' B \0, 0; t') ~ rig. This is because far 
into region I such locations r' B are beyond diffusional range of the interface: hence density 
correlations at r' B cannot have been influenced by reaction events during (0, t'). On the 
other hand, in region II (x' B < Xf), this conditional density will be strongly influenced by 
such reaction events. Whatever this density field may be, we expect that its maximum will 
never be greater than a value of the order of 71%°. Reactions tend to reduce densities, but we 
do not exclude the possibility that subtle B-A-B correlations could locally elevate the field 
somewhat. 

Let us now translate the above general physical expectations into two specific assumptions 
on the conditional density field. Simultaneously we introduce a third assumption, concerning 

P*A B {t). 

Assumption 1. There exists a positive finite constant U, such that: 

PbabWOAQ <v {20) 
n B 

This amounts to assuming that irrespective of what reaction-induced correlations exist be- 
tween points and r' B , the conditional density of B particles at r' B will always be less than, 
or at most of the order of, the far-field density of B reactants in the B bulk. 

Assumption 2. There exists a positive finite constant L, such that: 

P^B^m >L fo r ^ > i . (21) 

n B x t i 
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Roughly speaking this amounts to assuming that points in region I are uncorrelated with 
the interface. 

Assumption 3. p s AB {t) * s a decreasing function of time which is asymptotically a power 
law. 

Assumptions 1 and 2 immediately imply the same two assumptions but with A and B 
interchanged (since A and B are arbitrarily chosen labels). That is, the field paba{t'a\®i 0; t') 
appearing in obeys two assumptions analogous to 1 and 2. 

Based on these assumptions, we show in Appendix B that the contribution to I^(t) from 
integration over region I is a fraction of order unity of the value of I^(t). ( Equivalent ly, 
a fraction of order unity of the A-B interface pairs which involve one previously reacted A 
and are subtracted off by XI^(t), involve a B member originating from region I.) Therefore, 
since assumptions 1 and 2 imply that pbab equals rig in region I to within a finite pref- 
actor bounded above and below, it follows that if we replace Pbab{y'b\®-> 0> — ¥ n< B m the 
integrand of I^(t) in eq. (19), the result will equal the actual value of I^(t) to within a 
(time-dependent) prefactor of order unity, a(t). Making this replacement, using eq. fll3|), 
and performing an analogous replacement for I^(t), one obtains 

= ait) n% fdt'S^\t - t')p AB (t') , I*(f) = (3{t) n% f dt> S {1 \t - t')p s AB (t') , 

Jo Jo 

(22) 

where a(t),/3(t), are bounded positive functions of order unity, 

"min <a(t)< a max , p min <(3(t)< /5 max . (23) 

Here a min , a max , m \-n, Anax are finite positive constants. The one- dimensional return proba- 
bility S^(t) is defined as 

S*\t) = [ <b$($\o,$) n ±- . (24) 
J x t 

It measures the probability a reactant initially at the interface returns to the interface 
after time t. The scaling form, ~ is easily derived from the scaling form of the 

propagator \ eq. (p]). 

4. Reaction Rate in Compact Case (d + 1 < z) 

Having expressed the many-body terms I^(t), I^(t) in terms of the interfacial reactive 
pair density p s AB (t), we can solve for the reaction rate per unit area 1Z t = Xp s AB (t). According 
to the results of the previous section, eq. (p2f), the self-consistent solution for p AB (t), eq. 
(|15D, can be written 

p AB (t) = nfn'g - A fdt'S^\t - t')p s AB (t') - Xn(t) f dt> S {1 \t - t')p s AB {t') , 

JO JO 

n{t) = a(t)n£ + P(t)n% . (25) 

This "solution" of course involves the unknown function n(t). From the arguments of the 
previous section following from our assumptions 1 and 2, we know that n(t) is bounded 
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above and below. Now according to assumption 3, asymptotically p s AB (t) ~ t~ 5 with 5 > 0. 
Substituting this power law in eq. (25), and substituting x t ~ t x l z in the scaling forms of 
S^ d+1 '(t) and S^\t) from eqs. ( |I6"D and (|2^), one finds that as t — > oo the many-body term 
dominates over the other time-dependent terms in eq. (25), and up to a constant prefactor 
is equal to n(t)it 1 ~ 1 / 2: ~ <5 . It follows that at long enough times the many-body term must 
equal the first term on the rhs of eq. (25), n^n^, plus higher order corrections. Since n(t) 
is bounded, this implies that 5 = 1 — 1/ z and that n(t) tends to a constant at long times, 
n(oo). We will prove in section 7 that this constant is none other than the reactant density 
in the more dense of the two phases: 

n(oo) = . (26) 



oo 

l A . 



We remind the reader of our convention throughout this study, > nl 

Laplace transforming eq. (25), t —* E, and recalling that 1Z t = -VabW^ ^ * s simple to 
obtain the following self-consistent relation for the Laplace transform of the reaction rate 
per unit area, lZ t (E): 

* (F) = = n(EMS^(E)n t (E)} 

t{ } Eil + \S^(E) + \n% 7 (E)SU(E)}' 71 ; n%SW(E)lit(E) ' K ' 

Here, * indicates convolution in Laplace space. The function j(E) has a simple form for 
small E; since, by virtue of eq. (EBJ) , n(E — > 0) = n^/E, then from eq. ( BTJ ) one has 

7 (E) = l + 0(£) . (28) 

Now from eqs. ([I]), (|TB|) and S'^ +1 ^(t) and S^(t) are algebraic in time. Their Laplace 
transforms have the form S^ d+ ^{E) ~ i?^ 1 )/ 2 - 1 (valid only in compact dimensions, d+1 < 
z) and S^'(E) ~ (always valid). This section is concerned with the compact case; 

then we can rewrite TZt(E) in two ways: 



n t {E) 



E{l + (Etlf^ 1 ^ 1 + 7 (£) {Et* m i l/z) - 1 ) 



Xn°?n 



oo 
A lb B 



e[i 



(Et* 2 ) [{d+l)/z] ~ l \l+j(E) (Et 



-d/z 



(29) 



where 



*' = *" l^J ' C = (gtanfa-J ' ^ = ta {QfaJ 

(30) 

are essentially the three naturally occurring timescales introduced in section 1, generalized 
to the case of unequal initial reactant densities (ng 5 > n A ). It is important to note that the 
characteristic density determining these timescales is that of the denser bulk phase B. In eq. 
(29) for simplicity we have neglected numerical prefactors in the terms in the denominator. 
Note that the three characteristic timescales obey 

t* m = (t* 2 ) 1 - d/iz - 1 \t l ) d/{z - 1) (31) 

which implies that the magnitude of t* m always lies between those of t% and ti. Hence there are 
only 2 cases (see fig. ||). (a) For strongly reactive ("strong") systems, Q > Q*, the ordering 
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of timescales is t\ < t* m < t h (b) For "weak" systems, Q < Q*, one has t t < t* m < t\. The 
boundary between strong and weak regimes is defined by a critical effective local reactivity, 
at which t\ = t* m = t t : 

Q*t a = {n^a d )^ d+1 ^ d . (32) 

(a) Strong Systems : Q > Q*, t\ < t* m < t\. Before evaluating 1Z t in different time 
regimes, we note that the many body term (the third term in the denominator in eq. (29)) 
is unimportant whenever Eti ^> 1 (corresponding to t <C £;). Consider the term (that 
involving S^) in eq. (25) from which this many body contribution is derived. Now imagine 
replacing n(t) in this term by its maximum value, n(t) — > n max = « m ax n s + Pmax n A, sucn 
that n(E) = n max / E and hence j(E) = n max /n(oo) « 1. It would then indeed follow that 
the many body term in eq. (29) is higher order for Et\ ^> 1. Clearly, then, this must always 
be true. 

Consider firstly short times, E^ 1 <C t\_. The many body term can then be neglected. 
Considering the two cases E^ 1 <C and t\ E^ 1 ti, respectively, inverse Laplace 
transformation of eq. (29) yields 

r a (t<^) 

U t = kV>n%n ( g , k i2) w <^ . (33) 

( dx d t +l /dt ~ t^l*- 1 {t* 2 < t < t t ) 

These are 2nd order rate kinetics. An initial MF regime is followed at t* 2 by a DC regime. 
Notice that the timescale t* m is irrelevant. We remind the reader that our analysis does not 
describe times less than th (see comments following eq. (12)); hence, if Q is so great that 
eq. (|30D implies t\ <t h , then eq. (^) correctly describes the second DC regime only. 

Now consider very long times E" 1 ^> tf, the many-body term in eq. (29) is then dominant. 
Since at long enough times we may replace j(E) — > 1 as discussed, one now finds first order 
kinetics: 

7tt = JfeWrtf , k{l) ~^tt~ tl/Z ~ l ■ ^ 

Thus, at long enough times the reaction rate is controlled by the diffusion to the interface 
of the more dilute A species. 

(b) Weak Systems: Q < Q*, t\ < t* m < t\. Now the many-body term in the curly 
brackets in eq. (29) is much smaller than 1 whenever E^ 1 t* m (this can be seen by 
replacing n(t) with its maximum value as we did for (a) above). But for such E values, it is 
automatically true that Et\ <C 1 by virtue of the definition of weak systems {t\ > t£j, and 
hence the 2-body term in eq. (29) is also much smaller than unity. It follows that MF 2nd 
order kinetics pertain for all times less than t* m 

Tit = kWn<%n% , fc (2) w A (t < t m ) . (35) 

Notice that the 2-body and many-body terms are both proportional to negative powers 
of E, and that the magnitude of the many-body term's exponent is the greatest of the two. 
Now consider E^ 1 ^> t* m , when the many body term is much bigger than unity. It follows 
that this term is then also much bigger than the 2-body term, because t* m < t\ for these 
weak cases. Thus kinetics are first-order for t^> t* m : 

Ttt = k^n% , fc« « ^ ~ (t » t*J . (36) 
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To obtain eq. fl36|) we have replaced *y(E) — ► 1 for small i£. For weak systems, neither i£ 
nor ti are relevant. The reactivity is so small that the two-body term is never relevant, and 
2nd order DC kinetics are absent. 

5. Reaction Rate in Marginal Case (d + 1 = z) 

The previous section dealt with low compact dimensions, for which the 2-body return 
probability in Laplace space had the form S^ d+1 '(E) ~ £ , ( d+1 )/ z - 1 . For d+ 1 > z, this is no 
longer true. In this section we consider the marginal case, d + 1 = z; thus S^ d+1 ' ~ l/xf +1 ~ 
1/t, giving 

» jT dte-™-^ « ln[l/^] (Et & « 1) . (37) 

We have introduced a cut-off at t = t^, at shorter times S ,( - d+1 - ) (t) crosses over to a form 
appropriate to a <i-dimensional bulk problem, S^ d+1 ^ ~ l/(hxf) whose time integral gives a 
contribution of the same order as that from the lower limit in eq. ( p7|) . 

Aside from this modification, all steps leading to eq. (29) of the compact case are un- 
changed: the expression for the reaction rate IZt(E) (eq. (|27|) ) remains valid, and the form 
of (E) is unchanged. Thus, 

ME)- ! b 



E{i + l ^^+i(E)(Et* m y/^} 



1 (3£ 



^ \ 1 + ln(eQ/t h ) L 1 + 71^ J ln(VBt fc ) J J 

Here the definitions of t\ and are unchanged from the compact case (eq. fl3U|)), but now 

tl = {t h /e)e 1 ^ . (39) 



Let us define T} to be the time such that for E < T t 1 the many-body term (oc S^'(E)) 
dominates over the two-body term (oc S^ d+1 \E)) in the denominator of eq. (38) : 

^ = [ln(eT,/t,)] z/d . (40) 
w 

(We have included factors of e in the definitions of t\ and TJ above simply to ensure continuity 
of reaction rates; see eqs. ( fi^D and ( ^3|) below.) 

Analogously to the compact case, the condition t\ = t* m = T\ defines a critical reactive 
strength Q*, 

Q*t a = 1 , (41) 
In [eTi/t h \ 

defining the boundary between "weak" and "strong" kinetics (see fig. ^), for which it can 
be shown that the 3 relevant timescales have the same orderings as for the compact case. 

(a) Strong Systems : Q > Q*, t* 2 < t* m < T t . Consider first short times, E' 1 < T h 
Similar reasoning as for the compact cases implies that the many-body term in eq. (38) can 
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be neglected for such E values. Considering the two cases E 1 <C t% and E 1 ^> one 
obtains 

r a (t<zt* 2 ) 

1l t = k®n%n'g , fc (2) w I (42) 

( h d+1 /[t h \n(et/t h )} (Q<t<r,) 

The logarithm arises after inverse Laplace transformation of l/{£'ln(l/£ , t/ l )} which gives 
l/ln(i/i/j) for t ^> t^. This is shown in appendix C. 

For long times, E^ 1 ^> Tj, the many-body term dominates. Using 7(i£) ~ 1 for small 
enough which is easily demonstrated using similar arguments to those for the compact 
case, one finds first-order DC kinetics which are no different in structure to those for the 
compact case (see eqs. ([34|) and fl36|)): 

li t = k^n%, fed) » ^ „ t V*-i (t»T,) . (43) 

(b) Weak systems : Q < Q*, T\ <t* m < t\. For small times, E" 1 ^ t* m , the many body 
term is much less than unity; this is also true of the 2-body term since t\ > t* m (definition 
of weak system). On the other hand, when E~ l 3> t* m , the many body term is much larger 
than unity; it is also much bigger than the logarithmic 2-body term since E^ 1 ^> TJ follows 
automatically, because T] <t* m . Thus 

k^n^n™ , « A {t < t*J 

Kt={ . (44) 



6. Reaction Rate in Noncompact Case (d + 1 > z) 

In this section high non-compact dimensions are considered, d + 1 > 2. Perhaps the 
commonest physical example of small molecules [z = 2, d = 3) belongs to this class. Mathe- 
matically, the only distinguishing feature is that the Laplace transform of the 2-body return 
probability is now dominated by small times, since S^ d+lS) (t) of eq. (16) now decays faster 
than 1/t for times t > t^. 

, , roo I /+ \ (d+l)/z f 

S (d+1 \E) « dte-*-^ ) « ^ (Et fc « 1, * < d+l< 2+1) . (45) 

The above result is determined by the dominant cut-off at t — th. In fact it is valid only 
provided d < z because only then is the t < th time integral dominated by its upper limit, 
t — th- for these smallest times (which have been neglected in the original statement of our 
model, eq. (12)) one has in effect an infinite bulk reaction problem. It is as if the interface 
were infinitely large. Correspondingly, the true return probability is S*^ +1 ^ ~ l/{hxf) for 
t < th. When time integrated, for dimensions so high that even bulk reaction kinetics are 
non-compact, d > z, the lower cut-off at t a is now dominant, J t h e~ Et S^ d+1 ^dt P3 t a S^ d+1 \t a ) 



for Eth <C 1. This contribution now exceeds that displayed in eq. ([45]), and one has 



S^ d+1 \E)^^ (Et h <^l,d>z) . (46) 
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Consider firstly z<d+l<z+l. The reaction rate in Laplace space of eq. Q27T) now 
reads 

Xr) 00 ^) 00 

} ~ E{l + Q b t a (a d t h /h d t a )+l(E)(Et*JV*-i} " 1 ) 

Here j(E) is the quantity defined in eq. ( G7p and, as for the compact and marginal cases, 
it can be shown that j(E) « 1 for small enough i£. Thus for <5&i a < (a/h) z ~ d 

( kWn%n% , A;^ ps A (t <C O 

7t* = J ' (48) 

whilst for Qbt a > (a/h) z ~ d one has 

f k^n^n^ , k® « (t < t h {n^h d ) z ^ 1 - z '>) 

7^t= ' " (49) 

[ fc«nf , « dx t /dt ~ tV*-i (t > t h {n^h d ) z l^-^) 

Now consider the highest dimensions, d > z. Then 7Z t (E) is as in eq. fl47|) , except one 
replaces Qbt a ( ad th/h d t a ) — >■ Qfet a . This leads to eq. ([48]) which is now valid for all Q b values. 

7. Density profile 

We have seen that short time 2nd order reaction kinetics cross over at a regime-dependent 
timescale to 1st order diffusion-controlled kinetics. This suggests that the density fields on 
either side of the interface, ua^a] t), n B (r B ; t), are uniform for shorter times but develop 
depletion holes at the interface of size x± when the 1st order kinetics onset. To demonstrate 
this explicitly, we begin by using Doi's formalism in appendix A to derive the density field 
dynamics for small molecules (z = 2). We will then generalize results to arbitrary dynamical 
exponent z. For z = 2, we find 

|^ -DV A ^n A {vA]t) = -X5(x A )p S AB(t) , |^ - D V B \n B (r B ;t) = -X5(x B )p s AB (t) . 

(50) 

The sink terms on the right hand sides of eq. ( pT|) are proportional to the number of 
A — B pairs which are in contact at the interface, per unit area. Noting that translational 
invariance parallel to the interface plane implies n A ,UB depend on xa and xb only, eq. ( 50" ) 
has solution 

n A (x A ;t) = n%-\ [' dl/Gi%(x A )p AB (t r ) , n B (x B ;t) = - \ f dtGf} t ,{x B )p s AB {t') , 
Jo Jo 

(51) 

where Gf } (x) = / dr T G ( t 1] (r, 0) is the weighting for a particle, initially at the interface, to 
be distant x from the interface after time t. For arbitrary z, one just uses the appropriate 
propagator G t in eq. (|5T|). 

Before proceeding, let us use the above dynamics to prove eq. (ITOf) , 1Z t = \p s AB (t), a 
result that we have so far assumed as physically obvious. Now the total number of reactions 



per unit area \slZ t = J dx A [n A —u a (xa', t)]; integrating the first of eqs. (plD over all xa and 
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using the fact that G t (x) is normalized to unity, one has Jl t = \^dt'p s AB (t') which proves 
the desired result. 

In the below we need calculate only one of the density fields, say the less dense field n A , 
since one field uniquely implies the other. This follows after subtracting the two equations 
in (|5TD, and using G[ 1 \x) = Gt(—x), giving 

n A (x;t)-n B (-x;t)=n^-n^ . (52) 

That is, the difference between mean A and B reactant densities at equal distances from 
the interface is constant in time. 

7.1. Long Time Density at Interface 



It may appear that eq. fl5T|) together with eqs. (|T0| ) and (p7| ) provide a closed solution 



for the density fields. However, in fact eq. (|27| ) involves the unknown function j(E), whose 
small E behavior is needed to obtain the long time density fields. Hitherto we have asserted 
that its asymptotic behavior is 7(0) = 1, equivalent to the assertion that n(oo) = (see 
eqs. (PE|), (p^j) and fl2"8| ) and surrounding discussions). We must now prove these assertions. 
To do so, we will first argue that the A density at the interface, n s A (t) = n A (0; t), vanishes 
for long times. This extra piece of information will allow the determination of 7(0). 

We are able to prove n s A (oo) = by first relating n s A to the like particle correlation 
functions, p AA (t) = p AA (0,0;t) and p BB (t) = p BB (0,0;t), on the strength of the following 
physically motivated assumption on these functions: 

Assumption 4. 

[n A (t)] 2 < p s AA (t) , [n%(t)Y < p% B (t) . (53) 



This states that reaction-induced correlations can only increase density-density correla- 
tions of like particles, relative to the totally random case where one would have p s AA {t) = 
[n s A (t)] 2 . That is, we admit the possibility of clustering of like particles. 

To obtain information about p s AA and p s BB we first relate them to p AB - In Appendix D 
we use Doi's framework to derive dynamics for p AA , Pbb from which we derive the following 
exact equation 

P AA (t) + P s BB (t) = (n%-n%) 2 + 2p AB (t) + 2\ f dt> S^ d+l \t - t>)p AB {t') . (54) 

Jo 

According to the results of sections 4,5 and 6, at long times p s AB {t) = 1Z t /\ ~ v^'*' 1 . (Note 
this conclusion followed from assumption 3 and is quite independent of the numerical value 
of 7(0).) Substituting this power law in eq. ([54]), using S^ d+lS) {t) ~ t~( d+1 ^ z from eq. (|16D 
and incorporating cut-offs in the marginal and non-compact cases (see eqs. (|37|) , (f45|) and 
(p6|) ) one sees that the time-dependent terms on the right hand side of eq. fl54] ) tend to zero 
at long times. Thus, making use of eq. fl52|), we obtain from eq. ([54]) 



[A(oo) - 1] [n s A {oo)} 2 + [B(oo) - 1] KM] 2 = -2n^(oo)n^(oo) , (55) 
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where we have defined the unknown functions A and B such that p\^{t) = A(t)[n s A (t)} 2 and 
p s BB (t) = B{t)[n B {t)] 2 . Note that assumptions 4 imply A(t), B(t) > 1; this in turn implies 
that the lhs of eq. fl55|) is positive or zero. But the rhs is negative or zero. It follows that 
both sides of this equation must vanish, i. e. either or both of n s A (oo) and n s B (oo) vanish. But 
from eq. (0), n A (oo) < n s B (oo). Hence n s A (oo) = is proved. 



7.2. Full Density Field 



Having determined that n A (oo) vanishes, we return to eqs. (|5T| ) from which we will first 
determine 7(0) and then calculate the full density profile. Using the expression for 1Z t (E) 
in eq. (|27|) , and making the substitution p s AB (E) = TZt(E)/X, eq. (|5*1|) can be written in 
Laplace space as 



n: 



n A (x A ;E) = 



1 



Xn 



B 



1 + XS( d+1 \E) + An£ i{E)SW{E) 



Here the Laplace transform of the propagator Gf (x) has the following structure: 



(56) 



G%\x) = S {1) {E)g{xE 1 




(57) 



where g is a scaling function with the stated limits. We have used eq. ([14]) and the fact 
(see eq. (||)) that S^{t) = Gf\x = 0). 

We can now prove j(E = 0) = 1. Consider the limit E — > of the expression in eq. ( |56|) 
evaluated at = 0. In this limit the square bracket must vanish since n A (t = 00) = 0. Now 
for small enough E, the many body term Xn B 2r y(E)S( 1 > (E) always dominates over the other 
two terms 1 and XS^ d+1 \E) (see eqs. (29), (38) and (g7p). Thus, using G ( B \o) = S^{E), 
we must have 7(0) = 1. 

Consider now general values of xa, t and let us compare the two terms in the brackets on 
the rhs of eq. (|56"D . According to eq. ( f>% , the numerator of the 2nd term is less than or equal 
to the many body term in the denominator; it follows that the quotient can be comparable 
to 1 only for E values sufficiently small that the many body term dominates. As we saw in 
eqs. (29), (38) and (^7[), this corresponds to times longer than the timescale signifying the 
crossover from second to first order kinetics. Therefore, retaining leading order terms only 
in eq. (p§), one has 



where 



\x t 



n A (x A \t) 




' "short" times) 



/(«) 



00 



(u > r 



(51 



(59) 



and C~ l denotes inverse Laplace transform. Here by "short" times, we refer to times when 
second-order kinetics are valid. This completes our calculation of the density profile, which 
evidently confirms the physical expectations. One sees that at short times riA{xA',t) retains 
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its equilibrium value, whereas at longer times a reactant density depletion hole of size x t 
develops at the interface (see fig. ^). 

As an example of the form of f(u), consider small molecules (z = 2) for which G^p(x) 
is a Gaussian. Determining Gg (xa), one finds from eq. fl59D that f(u) = Erf(u); this 
is identical to the asymptotic density profile in the situation in which initially uniformly 



distributed small molecules adsorb irreversibly onto a surface [17,18 



8. Segregation Effects and Decay of Interfacial Density, n s A (t). 

We found in Section 7 that the reactant density at the interface, n s A (t), tends to zero at 
long times. In this section we determine the long time power law decay of n s A (t), considering 
for simplicity the symmetric case only, = . Interestingly, we will find segregation of 
reactants into A-rich and B-rich domains at the interface for the compact case. 

We begin by establishing time-dependent bounds on n s A {t). Now assumption 4 will lead to 
an upper bound of this type, because for the symmetric case p s AA (t) can be determined from 
eq. fl54|) since p s AB (t) is already known. What we need, in addition, is a lower time-dependent 
bound, which we now introduce by making one further assumption. This assumption is 
motivated by the physical expectation that the density of A-B interfacial pairs, p s AB {t), will 
never exceed the value it would have if there were no correlations between A and B particles, 
namely [n^(t)] 2 . That is, A-B reactions will always tend to diminish this pair density relative 
to the uncorrelated value. 

Assumption 5. There exists a positive finite constant b, such that: 

b P AB(t)<[n A (t)] 2 (n% = n%) . (60) 



(a) Noncompact Case (d + 1 > z). In Appendix D, eq. flD4|), using eq. ( p4j ) we show 



that for large enough times p s AA (t) ~ (n A /\)dx t /dt. Meanwhile, eqs. ( pRf ) and (^) imply 
that p s AB (t) ~ (n A I \)dx t /dt. Thus the upper and lower long time bounds on n s A implied, 
respectively, by assumptions 4 and 5 are proportional to the same algebraically decaying 
function of time. Hence 



nm*^^r~t^*»M (f-oo, d+l>z), (61) 
up to a (time-dependent) pref actor of order unity. 

(b) Compact Case (d+ 1 < z). For this shown in appendix D, eq. ([54]) leads 

to the conclusion (see eq. (|D5| )) that asymptotically p s AA (t) ~ n A x^ d ~ t~ d / z . Meanwhile, 
eqs. (|34l) and (|36"D imply the same decay for p s AB as for the noncompact case, p s AB {t) w 
(n%/\)dx t /dt ~ tC 1 "*)/* Hence the upper and lower bounds on n s A (t) implied by assumptions 
4 and 5 involve different power laws: n s A {t) decays at least as slowly as t~ d /( 2z ) and at least 
as rapidly as t( 1-z )/( 22 ). This is insufficient to determine the actual decay. 

We can make progress, however, by invoking the interface analogue of the arguments 



which were used by Ovchinnikov and Zeldovich |27| , and Toussaint and Wilczek |28[ to 
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analyze the bulk reaction system A + B — > 0. According to this generalization, which we 
have presented in the introduction, the density of A reactants at the interface cannot decay 
faster than \Jn c A Xt d ^ 2 ~ t~ d ^ 2z \ which is the rate determined by the decay of fluctuations 
in the random initial reactant distribution. But we have already shown that t~ d /( 2z ) is an 
upper bound. Hence 

n s A (t) « sjrf£xt d/2 ~ t~ d/(2z) (t -> oo , d + 1< z) . (62) 

Therefore, the asymptotic density decay at the interface is controlled by the rate of decay 
of fluctuations. It follows that A- rich and 5-rich regions of linear size Xt develop adjacent 
to the interface. These are illustrated schematically in fig. [7]. An important point to stress 
is that the long time reaction rate is itself not influenced by this segregation, to leading 
order: the long time reaction rate is governed merely by the fact that n s A (oo) = 0, whilst 
segregation effects are associated with higher order terms in n s A (t), i.e. the manner in which 
n s A decays to zero. 

9. Discussion 

We have shown here that the critical dimension for reaction kinetics at a fixed interface 
is d c = z — 1 . This is quite different to the result for reactions at a movable and broadening 
interface separating 2 miscible phases, which problem has been widely studied for the case 
z = 2 where d c = 2 has been found [|],^]|||9|]. For the fixed interface problem studied here, 
one has instead d c = 1. The difference between these 2 critical dimensions is due to the fact 
that for the case of miscible reactants, A-B reactions are not restricted to occur only in a 
(d — l)-dimensional plane. 

The most novel feature to have emerged from this study is that interfacial reaction kinetics 
are not of fixed order. This is rather unusual. For example, trimolecular, bimolecular and 
unimolecular reaction processes are generally governed by 3rd, 2nd and 1st order kinetics, 
respectively. The peculiar feature here is that 2nd order reaction rate laws are obeyed at 
short times, whilst 1st order kinetics describe the long time behavior: 

U t = k^n^n^ (2nd order) , K t = k {l) n^ (1st order) . (63) 

The 2nd order coefficient k^ may either be a constant (mean field kinetics, MF) or 
time-dependent (diffusion-controlled kinetics, DC). The time-dependence in the latter case 
is k^ f« dx d+1 /dt. In contrast, the 1st order kinetics are always DC, and the 1st order 
coefficient k^ 1 ' ~ dxt/dt is always time-dependent. 

An important feature of these 1st order kinetics concerns the different roles played by 
the two far- field reactant densities, n A and , in the case where they are unequal. The 
timescale at which these kinetics onset (either t* m or t{) is determined by the greatest, n# . 
However, the rate law itself involves the smallest one, n A (see eq. (|63|)). Correspondingly, 
in the region within a distance Xt of the interface the density profile falls to a value close to 
zero on the dilute A side, whereas on the denser B side the profile in this region drops to a 
finite value close to — n A . 

Apart from our main concern, the reaction rate, this paper has also addressed the evo- 
lution of density fields, key features of which are the densities at the interface n A , n s B . This 
enabled us to examine the validity of the "local mean field" decoupling approximation, 
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lZ t ~ Xn s B (t) n s A (t) , which approximates the densities on either side of the interface to be 
independent of one another, p s AB n s A n s B . In this picture the denser B side is viewed as 
presenting a uniform "reactive surface" of strength Xn s B (t) to the more dilute A reactants. 
Consider long times, when the diffusive flux of A particles per unit area at the interface is 
n A dxt/dt. Equating this to the reaction rate, one sees that the "local mean field" approxi- 
mation suggests 

n A^) ~ ^ s ~~df n A (t — > oo , local mean field approx.) (64) 

Now in the symmetric case, n s A — n%, the above result implies n s A {t) ~ t( 1_2 )/( 22 ). But we 
saw in section 8 that the decay rate is always limited by the rate at which fluctuations in the 
initial differences between densities on the A and B side near the interface can diffuse away. 
This limiting decay was shown to be ~ £~ d /( 2z ). This suggests that only in high dimensions, 
d+1 > z, is the n s A (t) ~ £( 1_2! )/( 22! ) prediction correct; indeed, we demonstrated this in section 
8. We conclude that the local mean field approximation is essentially valid for d + 1 > z 
at very long times. For all low dimensions d + 1 < z, however, fluctuations determine the 
decay law: n A (t) ~ t~ d '^ z \ segregation occurs at the interface, n s A and n s B are no longer 
independent and the the local mean-field approximation is wrong. 

Let us make a few comments about the interfacial densities in the asymmetric case, n s B > 
n s A . In this case we expect eq. ( |5"4] ) to be valid for all dimensions, since n^(oo) = — n A 
is then non- vanishing (see eq. (^)). Hence the B side will indeed supply a uniform reactive 
surface for the A reactants. Thus, we expect a different decay law, n s A ~ t( 1_2 )/ z for all 
dimensions. If the initial reactant densities n^ng 1 are almost but not quite equal to one 
another, we expect the symmetric case results will be valid up to a cross-over time at which 
n s B {t) drops to a value close to its asymptotic value, rf§ — n A . Thereafter, eq. (|64D will 
correctly describe n s A . 

We stress that this study has concerned irreversible reactions. Thus an equilibrium state 
is never attained. The final state will be governed by saturation effects at the interface, 
which have not been considered here. As t — > oo, in principle a final state will be attained in 
which reaction product fills every available surface site. (In practice, however, the timescale 
for this state to be reached may be experimentally inaccessible POjn 



To conclude, consider a few specific examples. An important parameter determining 
the class of reaction kinetics is the dimensionless local reactivity, Qbt a - Perhaps the most 
useful relation to help one estimate its value for a given system is Qbt a ~ k hnlk / k^!^ where 
/c bulk ~ Qbd 3 is the bulk rate constant, i. e. the 3-dimensional rate constant which would 
describe A-B reaction kinetics if the molecules could react anywhere within the bulk (see 
introduction). Here k^^ ~ a? /t a ~ 10 9 (litres/mol)sec _1 is the same quantity for radicals 
which are nature's most reactive chemical species. We assume here the molecular size a is 
roughly the same (a ~ 3A) in all small molecule cases. Thus if one has access to fc bulk for an 
A-B system, then one can estimate Qbt a - In the case where the reactive groups are attached 
to polymer chains, k hulk refers of course to the small molecule bulk analogue reaction system, 
i. e. the rate constant describing reactions between the same species after removal from their 
host polymer chains. 
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Small molecules: z — 2, d — 3. 

Consider firstly unequal initial bulk densities, n A ^ n^. The early 2nd order behavior 
is non-compact (d + 1 > z) and MF 2nd order kinetics pertain with k^ = h(Qi,a 3 ). These 
continue until t* m = D/[h(Qba 3 )n B °} 2 , when first order kinetics onset with time- dependent 
rate constant k^ = D/ (Dt) 1 ^ 2 where D is the molecular diffusivity. Note that the cross-over 
time t* m is determined by the greater of the two far-field densities, . 

The density profile on the less dense A side is (to leading order) identical to the "reactive 
surface" situation, having a depletion hole of size (Dt) 1 ^ 2 . The A density at the interface 
decays for long times to zero as n s A ~ l/t 1 ' 2 . There is no "hole" on the more dense B side, 
though the density is reduced from its initial value over a region extending (Dt) 1 ^ 2 into the 
bulk and has the long time value — n A at the interface. The symmetric case, = ng 5 , 
is different: there are long time holes on both sides and the interfacial density decay is 
n s A ~ 1/t 1 / 4 . 

Typical numerical values are a ~ h ph 3A and D ph 10~ 5 cm 2 /sec. Now for the vast 
majority of reacting species, k hnlk ~ 10 3 (litres/mol)sec _1 , implying t* m ~ 10 sec/0^, where 
<Pb — n^a 3 is the far- field volume fraction of B reactants. Thus, depending on the value 
of 4>b, this timescale may become so large that the diffusion-controlled kinetics get washed 
out by other effects such as convection. For highly reactive species such as radicals, on the 
other hand, one has k hnlk ph 10 9 (litres/mol)sec _1 and t* m w 10~ 10 sec/</>^; these kinetics are 
then observable over a very large range of densities. 



Small molecules in d — 1. 

This is a marginal situation (z = d + 1) arising in systems where small molecules (z = 2) 
are restricted to an effectively one-dimensional geometry, e. g. molecules trapped in a thin 
tube. 

For highly reactive species, Qbt a ~ 1 (i- e. fc bulk ph k]?^) the initial regime is 2nd order with 
a weakly time-dependent rate constant k^ ph D/ \n(t/t h ). At time T x = D~ 1 (ng > )" 2 [ln(ng > /i) 2 ] 2 
first order kinetics onset with time-dependent k^ = D / \Dt) 1 / 2 . 

For most cases, however, the local dimensionless reactivity Qbt a will be below a very high 
threshold value (i.e. very close to unity) given by Qlt a = (a / h) / [\ia(n c g h) 2 ] 2 . In such cases 
an initial 2nd order mean field regime with k^ 2 > = Qbha = k hulk h/a 2 is followed at time 
t* m = ta/{Qbta n< Bh) 2 by the same 1st order kinetics. In all cases a depletion hole of size 
(Dt) 1 / 2 grows at long times on the dilute A side. 



Unentangled polymers, short times: z — 4, d — 3. 

Consider an interface separating two immiscible unentangled polymer melts comprising 
chains with degree of polymerization iV and radius of gyration R, each carrying one reactive 
group. Thus the density of reactive groups is n A = = l/(Na 3 ) or equivalently 4>a = 
cpB — 1/iV. Then Rouse dynamics apply [p5,p6|], x t ~ R^/t) 1 ^ for times less than the single 
chain longest relaxation time, r ph t a N 2 . Thus, for t < r, we are in the marginal situation 
z = d + 1 = 4. 

Consider first the maximally reactive case Qbt a = 1- Initially kinetics are 2nd order with 
weakly time-dependent rate constant k^ = (i? 4 /r)/ln(t/t/ l ). The cross over to first order 
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kinetics, with ~ (R/t) (t/r) 3 ^ 4 , occurs at TJ = t a ^ s ln(0 B 4 ' 3 a 4 //i 4 ). For typical 
values ft/a = 5, N = 200 one has 7} « 0.02 r. 

For less reactive species, Q&t a < Qj^a where ~ a/[^ln(0 B 4//3 a 4 //i 4 )], the kinetics are 
different. (For the above typical numerical values, Qlt a ~ 0.7.) In this case second order MF 
kinetics with k^ = h(Qba 3 ) are followed at t* m = t a [Qb(h/a)t a ^)B\~^ 3 by 1st order kinetics 
with jfeW = {R/T){t/r)- 3 / A . 

In both of these examples, a long time depletion hole grows on the dilute A side whose 
size increases in time as x t ~ R^t/r) 1 ^ 4 . 



Entangled polymers, "breathing" modes: z = 8, d = 3. 

Consider the same polymer example as above, but now chains are entangled. Using the 
reptation model to describe the polymer dynamics, let us ask what reaction kinetics are 



during the short time "breathing modes" regime (t e < t < £&) when fl25|f26|| x t = r e {t/t e ) ' . 
Here t e = N 2 t a is the entanglement time (N e being the entanglement threshold), tf, = 
(N/N e ) 2 t e is the Rouse time for the one-dimensional tube motion and r e = N^ 2 a is the tube 
diameter. This is an interesting example of a compact case, d + 1 < z — 8. 

Consider very reactive groups such as radicals, Qbt a ~ 1, and n# values such that t\ (the 
diffusion time corresponding to a distance equal to the typical separation between the B 
reactive groups) satisfies t e < t\ < Then ti = t^n^r 3 )" 8 ^ 3 . Now for t > t e there is no 
MF regime and kinetics are 2nd order DC, k^ « (r i /t e )(t/t e )~ 1 ^ 2 . For t > U the kinetics 
become 1st order with k^ = {r e /t e ){t/t e )~ 7 l 8 and a depletion hole grows on the A side of 
size ~ t 1/8 . 

As a specific example, if all B-chains carry one reactive end-group (n^a 3 = 1/N) and 
iV « 10 4 , N e « 200 one has t t « 30t e « 10~ 3 t 6 . Thus both 2nd and 1st order kinetics as 
described above will occur within the t 1 ^ 8 regime. 
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Appendix A. Derivation for z = 2 of Dynamical Equations for 2-Body and 1- 
Body Density Correlation Functions. 



In this Appendix we employ the second-quantization formalism for classical many-particle 



systems developed by Doi j42] to derive exact evolution equations, in the case of small 
molecules (z — 2), obeyed by the 2-body correlation functions p A B, Paa, Pbb and the density 
fields riA, n B . We do not attempt here a self-contained discussion of the Doi formalism: the 
reader is referred to Doi's papers, refs. for the necessary background. 

In the Doi formalism any physical quantity A is mapped onto a quantum operator A 
given in terms of the Bose creation and annihilation operators. For systems consisting of 
two types of particles, A and B, there are two types of Bose operators, ^(r), -^(r), for 
every spatial location r. These satisfy the following commutation relations 



[^(r),^t(r')]=5(r-r>, 



V[i 1 



t(r),4(r')] = [^(r) l ^(r0] = 



(Al) 



where v — A, B and fi = A,B. The dynamics of A(t) are determined by a quantum 
propagator G, 



A(t) 



Ae 



-Gt 



(A2) 



Here, (1| 



exp[J drA^AiyA) J dr B ip B (r B )} is a coherent state, (0| is the vacuum state 
= (1| . The quantum state |c) represents the initial state of the 



and (1| ip A (r 

system; although its form will not be relevant in the subsequent calculations, we remark that 
in the case where A and B particles are initially randomly distributed with densities and 

ng> then |c) = expfn^ng 1 / dr A ip\(r A ) J dr B ip B (r B )) |0). 

In the present reaction-diffusion interface problem, the propagator G = G Q + G r consists 
of a "diffusion" part, Gq, and a "reaction" part, G r . For small molecules obeying simple 
Fickian diffusion, according to Doi 



G = -D dr A ^ A (r A )V A M r A) ~ D / dr B ^ B (r B )V B ^ B (r B ) 



(A3) 



The reaction part, G T , is constructed from the reaction sink function which in our model 
is XS(x A )S(r A — r B ). We find that the corresponding quantum operator is 



G r = X I dv 1 ' { ^ A (r T )4(r T )Mr T )Mr T ) ~ ^(r T )^(r T ) 



(A4) 



The two annihilation operators in the second term in G r are responsible for reactions between 
A-B pairs in contact at the interface, while the first part is necessary to ensure proper 
normalization of averages ||42|. 

Differentiating eq. (|A^), and using the identity (1| G = which follows from the above 
definition of G, one has 



dA(t) 
dt 



[G, A}e 



-Gt 



(A5) 



We can now derive the dynamical equations obeyed by the correlation functions from eq. 
). The operator representations of the many-body correlation functions are [021 as follows: 



IK 
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/V(r,r') = vJ(r)V>I(rOVv(r)Vv( r ') > 
p^(r,r',r") = 4(r)^(r')4(r")^(r)^(r')^(r") (<t,/x,i/ = A, 5) . 

(A6) 



From eq. ( |A5| ) , with A = ua and A = tia, and using the above representation for ua, one 
obtains 

dnA{ *f t] = DV\ (l\A(r A )MrA)e- dt \c)-XS(x A ) <1| 4(^)4(r5)^(^)^(^)e- 5 ' |c> 

(A7) 

after using the commutation relations of eq. ( |A1[ ) and the properties of the coherent state 
(1|. From eqs. (A6) and eq. ( |A2| ) one recognizes this as the first of eqs. ( |50D in the 
main text (dynamics of ua,tib)- Notice that the commutation of G with n in eq. ( |A5| ) 
produced, among other quantities, a higher order correlation function; this is the origin of 
the hierarchical structure of the reaction-diffusion equations. 

Performing a similar analysis for the 2-body correlation functions (setting A = p^ v and 
A = p^u) one obtains eqs. (12) (dynamics of pab) and eq. (Dl) (dynamics of Paa,Pbb)- 

Appendix B. Proof That Relative Contribution of Region I to J A (t) is Order 
Unity 

We saw in section 3 that the many-body integral term I^(t) involving ps J 4s( r |0, 0; t) in 
eq. (19) receives contributions from two space-time regions, I and II. In this appendix we 
demonstrate that the contribution from region I is a fraction of order unity of the value of 
Z^(t) itself. We show this by firstly deriving an upper bound on I^(t). Then we derive a 
lower bound on J A (t) corresponding to Pbab being zero in region II and having its minimum 
value in region I. The lower and upper bounds will then be shown to be of the same order, 
proving the desired result. These bounds are all consequences of assumptions 1, 2 and 3 (see 
main text). 

From assumption 1 one sees that the substitution Pbab(?'b\®i 0; ~^ Un.% in the expres- 
sion for Z^j(t) (eq. (19)) defines an upper bound on I^(t), 

Ii{t) < Un% fdt'SV{t-t') P s AB {t') , (Bl) 
Jo 

where we have used eq. (|T^) to perform the integration. 
Now assumption 2 implies that 

Ln^ (x' B /x t i > 1 , region I) 
PBAB(r' B \0,0;t')>{ ' ' . (B2) 

(x' B /x t > < 1 , region II) 

Substituting the above lower bound on pbab in the expression for /„(£) of eq. (19), we 
obtain the following lower bound on I^(t) 

Ii{t) > Ln™ f dt> — p AB (t') f d d i h(0 , (B3) 
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where we have used the following scaling structure for the function F t (r' B ) appearing in the 
expression for I^{t) of eq. (19), 



Ft(r' B ) 



x t 



V x t , 



h(u) 



i {u < 1; 
o (n>r 



(B4) 



This follows from eq. (|14|). The integration variable in eq. ( |B3| ) is £ = r' B /xt-t' and £ x 
denotes the component of £ orthogonal to the interface. 

Now using eq. (|B4]), there exists a positive constant E of order unity such that 



E (xt/xt-t < 1 or H < t/2) 
d d Zh(£)>{ (B5) 
£,X«v/»t-f ) | o {x t ,j Xt _ tl > 1 or H > t/2) 



Expression 



) in inequality ( |B3| ) implies that 

t/2 



i-t/2 1 

J^(t) > £Ln£ / eft* p^(t') . 



(B6) 



Notice that according to eq. (0), S^(t) ~ Therefore, inequality ( [B6D is very close to 
showing that the lower bound on l£(t) is of the same order as the upper bound in eq. (|B1|) , 
except for the fact that the time integral on the right hand side of this inequality has upper 
limit t/2 rather than t. However, if one makes the replacement t/2 — > t for this upper limit, 
this yields the same result to within a constant of order unity. This is a consequence firstly 
of the fact that since z > 1 in eq. ([j]) thus both J ^ 2 dt' / x t -t' an d J^ 2 dt'/ x t-t' are of the 
same order, and secondly that p^b(£')) according to assumption 3, is a decreasing function 
of time. 

Thus we have shown that the upper and lower bounds on I^(t) are of the same order. 
But the crucial point is that the lower bound on I^it) of eq. ( p3|) results from an integration 
receiving zero contribution from region II and minimal contribution from region I. It follows 
that the actual contribution from region I must be of the same order as the actual value of 



Appendix C. Inverse Laplace Transform of l/(E\nE) 
We show in this Appendix that 



C 



-i 



-1 



E \n(Et h 



1 



r(l)ln(t/4) 



(t > t h ) , 



(CI) 



where £ 1 denotes inverse Laplace transform. Now it is well known that 



c- 1 



1 



(Et h ) 



in— 1 



tlT{n) 



[n > 0) . 



(C2) 



Integrating both sides of eq. ( |C2|) with respect to n from n = to n = 1 then gives 



c- 1 



1 






— 1 




K 1 Et h j 


In Et h . 


t Jo 



-dn 



(C3) 
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The function 1/T(n) is finite for all n in [0, 1]. Thus, expanding 1/T(n) in a Taylor series 
around n — 1 in the integral on the rhs of eq. (|C5|), we obtain 



1 e nln(t/t h ) 



t Jo T(n) T(l)t h \n(t/t h ) \t h [\n(t/t h )\ 2 J 

Considering the limit t ^> t h , corresponding to Et h <C 1, from eqs. (|C3|) and ( |C4]) we deduce 



eq. <&■ 



Appendix D. Like Particle Correlation Functions Paa-iPbb- Asymptotic Decay 
and Proof of Eq. (55) 



Using Doi's second-quantization formalism |42| it is shown in Appendix A that for z = 2 
(small molecules) Paa{ y a, r ' A , t) and ps^rg, r' B ; t) obey the following dynamics 

I- 

I dt 



al ~ D t V L + V r'J \ Paa(*A, r' A ; t) = -X5(x A )p AAB (r a , r' A , r A ; t) - X5(x' A )p AAB (r Al r A , r' A 



d_ 
dt 



- D [Vj B + Vr'J \ Pbb^b, r' B ; t) = -X5(x B )p BBA (r B , r' B , r B ; t) - X5(x' B )p BBA (r B , r' B , r\ 



Eqs. (Dl) have a similar form to eq. (12) for the pab{?a, ?b', t) dynamics except that the two- 
body sink term responsible for pair reactions in eq. (12) is absent since two particles of the 
same species cannot react with one another. Solving eqs. (Dl) and setting r A , r' A , r B ,r' B = 
we obtain 

P S AA(t) = «T ~ 2A/*(t) , p s BB (t) = {n™Y - 2XI^(t) , (D2) 

where I^(t), I^(t) are defined in eq. (17) of section 2 and are the identical many-body 
integral expressions which appeared in the solution for p s AB of eq. fll5|). 

Using the above expressions for the like interfacial pair densities together with eq. fll5D 
for p AB , one obtains the exact relation displayed in eq. fl54]) . This result tells us we can 
determine the sum p s AA (t) + p BB (t) from knowledge of p AB {t). For general z, eq. (|54]) remains 
valid, provided one replaces everywhere the Gaussian z = 2 propagator with the propagator 
G t appropriate to the dynamics. 

Consider the symmetric case, p s AA {t) = p s BB (t) • Laplace transforming eq. ([7J) we obtain 

p AA {E) = p AB (E) {l + XS^(E)} « = (D3) 

In the non-compact case (d + 1 > z), from eqs. (fl5|) and fl46|) , according to which S^ d+1 '(E) 
is a constant, and from the asymptotic form of p AB (t) = lZ t /X ~ t 1 / 2 " 1 (see eqs. ([48]) and 
(H)) one has 

f {(1 + \t h /h d+l )/X] n^dxt/dt (z<d+l<z + l) 
p AA (t) « « = n£, t -> oo) 

{ {(1 + Q b t a )/X} dx t /dt (d > z) 

(D4) 

Meanwhile in the compact case, from eq. ( |i~6|) one has S^ d+1 \E) ~ E^ d+1 ^ z 1 ; thus, using 
p s AB (E) = lit(E)/\ ~ E- X l z from eq. (29), valid for E -> 0, one has from eq. (|D3p 

« "a ^ ~ t' d/z {d<z , = n^, t -> oo) . (D5) 
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Figures 



FIG. 1. Two bulk phases A and B, separated by a thin interface of width h, contain dif- 
fusing reactants A and B with densities n^ng 1 . Reactions between A and B molecules, 
generating inert products, may occur in the interfacial region only. The typical distance 
between reactants on the B side is I = a{n B °a d )~ l l d . 

During short time 2nd order diffusion-controlled kinetics regimes, the reaction rate is 
determined by the small fraction of A-B pairs which were initially close enough to have 
diffused and met within time t. That is, reactions are confined to those pairs whose 
exploration volumes (indicated by dashed lines) overlap at time t. Note that such pairs 
must be within xt of the interface. 

FIG. 2. Schematic of the trajectory of an A particle after time t, given this particle was 
initially within diffusive range of the interface. Since the number of encounters with the 
interface is an increasing function of time, even for relatively weakly reactive species the 
A particle is certain to have reacted at sufficiently long times. The timescale is either t* m 
or ti (see main text). 

FIG. 3. A reactant density depletion hole of size Xt ~ t l l z grows at the interface for long 
times, (a) The symmetric case, n°£ = n^. The reactant density at the interface, n s A , tends 
asymptotically to zero, (b) In the asymmetric case, > n A , the interfacial density n s A 
on the dilute side tends to zero, while n s B asymptotes nf — n A . 

FIG. 4. The depletion in the number density of reactive A-B pairs at the origin from the 
value it would have in the absence of reactions originates from three terms, eq. ([H|). The 
two-body term counts those A-B pairs which would have been at the origin at time t, but 
failed to arrive because both members reacted at an earlier time t' at point . The first 
many-body term (l£(t), described in the figure) counts A-B pairs which would have been 
at the origin had there been no reactions, but failed to arrive because one member of the 
pair, the A member, reacted at an earlier time. The second many body term, I^t), is 
identical except the roles of A and B are interchanged. 

FIG. 5. Trajectories of a typical A-B pair which at time t is at the origin, and whose A 
member reacted at the interface at time t' with another B-type particle (B). Trajectories 
are shown projected onto the x — r plane where x is distance from the interface and r 
is time. The properties of trajectories of this type determine the value of the integration 
which determines the many body term I^it) (see eq. (19) and following discussion in 
main text). At the B particle of the pair has x-coordinate x' B which can be (a) in 
region I, if x' B > x T (as shown in the figure), or (b) in region II, if x' B < x T . 

FIG. 6. Reaction rate per unit area as a function of time in the Q-n^ plane (n^ > n^). 
Units are chosen such that a = t a = 1. (a) Compact case (d + 1 < z). (b) Marginal case 
(d + 1 = z). Reaction kinetics in the noncompact case are the same as in the "weak" 
regions of (a) and (b). 

FIG. 7. Schematic representation of asymptotic segregation of reactants into A-rich and B- 
rich domains of size x t near to the interface. This segregation occurs only for sufficiently 
low dimensions, d + 1 < z. 
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Fig. 1 
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Fig. 
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Fig. 3 
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Fig. 4 
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Fig. 5 



(a) (d+l)<z 
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Fig. 6(a) 



(b) (d+l) = z 
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Fig. 6(b) 
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Fig. 



